Non-integrability of the generalised 
spring-pendulum problem 



Andrzej J. Maciejewskif Maria Przybylska|§ and 
Jacques- Arthur Weil$^| 

f Institute of Astronomy, University of Zielona Gora, Podgorna 50, 65-246 Zielona 
Gora, Poland 

J INRIA - Projet Cafe, 2004, Route des Lucioles, B.P. 93, 06902 Sophia Antipolis 
Cedex, France 

§ Torun Centre for Astronomy, Nicholaus Copernicus University, Gagarina 11, 87-100 
Torun, Poland 

If LACO, Faculte de Sciences, 123 avenue Albert Thomas, 87060 Limoges Cedex, 
France 

Abstract. We investigate a generalisation of the three dimensional spring-pendulum 
system. The problem depends on two real parameters (k,a), where k is the Young 
modulus of the spring and a describes the nonlinearity of elastic forces. We show 
that this system is not integrable when k ^ —a. We carefully investigated the case 
k = —a when the necessary condition for integrability given by the Morales-Ramis 
theory is satisfied. We discuss an application of the higher order variational equations 
for proving the non-integrability in this case. 



PACS numbers: 05.45.-a,02.30.Hq,45.20Jj 

Submitted to: J. Phys. A: Math. Gen. 

E-mail: maciejkaOastro . ia.uz . zgora. pi, Maria. PrzybylskaOsophia. inria.fr, 
Jacques-Arthur . Weil@inria.fr 



Non-integrability of the generalised spring-pendulum problem 



2 



1. Introduction 

The spring-pendulum, which is known also under other names: swinging spring or 
elastic pendulum, is a very simple mechanical system having a very complex dynamical 
behaviour and this is why sometimes it is included to nonlinear paradigms. It consists 
of a point with mass m suspended from a fixed point by a light spring, moving under 
a constant vertical gravitation field. In Cartesian coordinates (x, y, z) with the origin 
at the point of suspension of the pendulum, the system is described by the following 
Hamiltonian 

H ° = + P ' + P ^ + m9Z + ~ Zo)2 ' (1) 

where r = x 2 + y 2 + z 2 , 1$ is the unstretched length of the spring, k G M + is the 
Young modulus of the spring. The motion of this system is a complicated combination 
of two motions: swinging like a pendulum and bouncing up and down like a spring. 

According to our knowledge, this system appeared first in [I] as a simple classical 
analogue for the quantum phenomenon of Fermi resonance in the infra-red spectrum of 
carbon dioxide. More about the history of this system can be found in j2j. Recently it 
has been analyzed in connection with the modelling of phenomena in the atmosphere 
[HI El Ej. Because of the complicated dynamics, various approaches for its analysis 
were applied: asymptotic methods 0, various perturbation methods [7J|E|, numerical 
methods 9j, various formulations of KAM theorem, the Poincare section, the Lapunov 
exponents JU]; the Melnikov method [TT) I12j. etc. A brief review of a large number of 
earlier papers on the spring-pendulum can be found in [2j and [T3j . 

Hamiltonian system generated by possesses two first integrals: Hamilton 
function H Q and the third component of the angular momentum 

p z = xy - yx, 

and for its complete integrabilty in the Liouville sense the third first integral is missing. 
Numerical computations suggest that such additional first integral does not exist and 
the system is chaotic. The first rigorous non-integrability proof for this system was 
obtained by Churchill et al |T3| by means of the Ziglin theory E>] . This result can 
be formulated in the following theorem. 

Theorem 1 // the Hamiltonian system given by Hamiltonian function Ho is integrable 
with meromorphic first integrals in the Liouville sense, then 



where q is a rational number. 

Morales and Ramis using their theory formulated in obtained a stronger result; 
they restricted the family (J2J) of values of parameter k for which the system can be 
integable. Namely, they proved in ^H] the following. 
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Theorem 2 // the Hamiltonian system given by Hamiltonian function H is integrable 
with meromorphic first integrals in the Liouville sense, then 

p(p + 1) 
p 2 + p — 2 ' 
where p is an integer. 



From the above theorem it easily follows that the physical spring-pendulum with 
k G M + is non-integrable except for the case k = 0. For k = the system is integrable 
because of separation of variables in the potential. 

In fact, the results presented above concern the two dimensional spring-pendulum 
system obtained in the following way. If we choose initial conditions in such a way that 
value of p z equals zero, then the motion takes place in a vertical plane and we obtain 
the two-dimensional system. But the non-integrability of the two-dimensional spring- 
pendulum implies immediately the non-integrability of the three-dimensional spring- 
pendulum. 

The aim of this paper is to investigate the integrability of the spring-pendulum 
system when the elastic potential contains also a cubic term. In other words, we consider 
a generalised spring-pendulum system described by the following Hamiltonian 

H = ^-(pl + p 2 y + p 2 z ) + mgz + -k{r - l ) 2 - -a(r - l ) 3 , (4) 

where k G M + and oGR. Our main result is the following. 

Theorem 3 // the Hamiltonian system given by Hamiltonian function (j3J) is integrable 
with meromorphic first integrals in the Liouville sense, then k = —a. 

In our proof of this theorem we apply the Morales-Ramis theory and some tools of 
differential algebra. Basic facts from the Morales-Ramis theory and some results 
concerning special linear differential equations are presented in Section 2. We derive 
variational equations and the normal variational equations for a family of particular 
solutions in Section 3. Theorem |3] is proved in Sections 4 (case a = 0) and 5 (case 
a / 0). In Section 4 we revise the result of Morales formulated in Theorem El Namely, 
we show that for values of k given by condition (JHJ) the system is non-integrable except 
for the case k = 0. In this section we also show two different kinds of arguments 
which give rise to non-integrability of the classical spring-pendulum system when k > 0. 
In Section 6, we study the exceptional case a = —k and conclude that the Morales- 
Ramis method yields no obstruction to integrability, whereas dynamical analysis seems 
to indicate that the system is not completely integrable. 

The Morales-Ramis theory was applied to study the integrability of many 
Hamiltonian systems, see examples in book [T2j and in papers [111 E01 EU EHJ EU EHl 
EEl E3 EE] ■ The differential Galois approach was used also for proving non-integrability 
of non-Hamiltonian systems, see [213 EH EH] • Difficulties in application of this theory 
can be of a different nature but mainly depend on dimensionality of the problem and 
the number of parameters. Although it seems that the Morales-Ramis theory gives the 
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strongest necessary conditions for the integrability, as far as we know, no new integrable 
system was found with the help of it. For the system investigated in this paper we have 
a very exceptional situation: we found a one parameter family of Hamiltonian systems 
for which the necessary conditions of integrability are satisfied, but, nevertheless, there 
is evidence that this family is not integrable, either. Another example of such family 
can be found in 

2. Theory 

Below we only mention basic notions and facts concerning the Morales-Ramis theory 
following [13132!. 

Let us consider a system of differential equations 

= v(x), teC, x e M, (5) 

defined on a complex n-dimensional manifold M. If <p(t) is a non-equilibrium solution 
of ©, then the maximal analytic continuation of (p(t) defines a Riemann surface F 
with i as a local coordinate. Together with system (J5J) we can also consider variational 
equations (VEs) restricted to XrM, i.e. 

£ = T(v)C, £ e T r M. (6) 

We can always reduce the order of this system by one considering the induced system 
on the normal bundle N := T^M/TT of V [33] 

i] = ir i! [T(v)7r~ 1 r)) , V eN. (7) 

Here n : TpM — > N is the projection. The system of s = n — 1 equations obtained in this 
way yields the so-called normal variational equations (NVEs). The monodromy group 
M. of system (JJJ) is the image of the fundamental group 7r 1 (r,t ) of T obtained in the 
process of continuation of local solutions of (J7J) defined in a neighbourhood of t along 
closed paths with the base point to- By definition it is obvious that M. C GL(s, C). A 
non-constant rational function f(z) of s variables z = (zi, . . . ,z s ) is called an integral 
(or invariant) of the monodromy group if f(g ■ z) = f(z) for all g e A4. 

In his two fundamental papers [HUE)], Ziglin showed that if system (J3J) possesses 
a meromorphic first integral, then the monodromy group M. of the normal variational 
equations ((Zj) has a rational integral (an invariant function). This result allowed him to 
formulate a necessary condition for the integrability of Hamiltonian system. 

If system (JHJ) is Hamiltonian then necessarily n = 2m and there exists function 
H on M such that u(v,u) = dH ■ u for an arbitrary vector field u on M (here u 
denotes a symplectic structure on M). For a given particular solution (p(t) we fix 
the energy level E = H(<p(t)). Restricting (jSJ) to this level, we obtain a well defined 
system on (n — 1) dimensional manifold with a known particular solution <p(t). For this 
restricted system we perform the reduction of order of variational equations. Thus, the 
normal variational equations for a Hamiltonian system with m degrees of freedom have 
dimension s = 2(m — 1) and their monodromy group is a subgroup of Sp(s, C). The 
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spectrum of an element of the monodromy group g G Ai C Sp(2(m — 1),C) has the 
form 

spectr(#) = (Ai, A^ 1 , . . . , A m _i, A" 1 ^), A, G C, 
and g is called resonant if 

m—1 m—1 

JjAf' = l for some (ki, . . . , fc m _i) G Z m ~\ ^ k, ^ 0. 
1=1 i=i 

In ^H] Ziglin proved the main theorem of his theory. Here we formulate it as in jHH] • 

Theorem 4 Let us assume that there exists a non-resonant element g G Ai . If 
the Hamiltonian system with m degrees of freedom has m meromorphic first integrals 
Fi = H, . . . , F m , which are functionally independent in a connected neigbourhood of 
T, then any other monodromy matrix g' G Ai transforms eigenvectors of g to its 
eigenvectors. 

Recently Morales-Ruiz and Ramis generalised the Ziglin approach replacing the 
monodromy group Ai by the differential Galois group Q of NVEs, see |17[ For 
a precise definition of the differential Galois group see [313 [T7J EE] • We can consider Q as 
a subgroup of GL(s, C) which acts on fundamental solutions of (|IJ) and does not change 
polynomial relations among them. In particular, this group maps one fundamental 
solution to other fundamental solutions. Moreover, it can be shown that Ai C Q 
and Q is an algebraic subgroup of GL(s,C). Thus, it is a union of disjoint connected 
components. One of them containing the identity is called the identity component of Q 
and is denoted by Q°. 

Morales-Ruiz and Ramis formulated a new criterion of the non-integrability for 
Hamiltonian systems in terms of the properties of Q° 01] . 

Theorem 5 Assume that a Hamiltonian system is meromorphically integrable in the 
Liouville sense in a neigbourhood of the analytic curve T. Then the identity component 
of the differential Galois group of NVEs associated with T is Abelian. 

We see that assumptions in the above theorem are stronger than in the Ziglin theorem. 
Moreover, as Q D Ai, Theorem 03 gives stronger necessary integrability conditions than 
the Ziglin criterion. 

In applications of the Morales-Ramis criterion the first step is to find a non- 
equilibrium particular solution, very often it lies on an invariant submanifold. Next, 
we calculate VEs and NVEs. In the last step we have to check if Q° of obtained NVEs 
is Abelian. Very often in applications we check only if Q° is solvable, because if it is 
not, then the system is not integrable. 

For some systems the necessary conditions for the integrability formulated in 
Theorem |S] are satisfied, but, nevertheless, they are non- integrable. In such cases, to 
prove the non-integrability we can use the stronger version of the Morales-Ramis theorem 
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based on higher orders variational equations [T71 137] . The idea of higher variational 
equations is following. For system (j2J) with a particular solution tp(t) we put 

x = <p(t) + + e 2 £ (2) + • • • + e k ^ k) + • • • , 

where e is a formal small parameter. Inserting the above expansion into equation © 
and comparing terms of the same order with respect to e we obtain the following chain 
of linear inhomogeneous equations 

d_ e W = A (t)£« + . . . , C^), fc = 1, 2, . . . , (8) 

where 

and /i = 0. For a given equation (JHJ) is called /c-th order variational equation (VE^). 
Notice that VEi coincides with There is an appropriate framework allowing to define 
the differential Galois group of k-th order variational equation, for details see [TT| I37|. 
The following theorem was announced in [37j . 

Theorem 6 Assume that a Hamiltonian system is meromorphically integrable in the 
Liouville sense in a neighbourhood of the analytic curve T. Then the identity components 
of the differential Galois group of the k-th order variational equations VE^ is Abelian 
for any k G N. 

There is also a possibility that the differential Galois groups of an arbitrary order 
variational equations are Abelian. Then we have to use another particular solution for 
the non-integrability proof. 

For further considerations we need some known facts about linear differential 
equations of special forms. At first we consider the Riemann P equation 
d 2 rj f 1 — a — a' 1 — 7 — 7 ; \ dr? 



(L: 2 V z z — 1 J dz ^ 

act 77' 13(3' — aa' — 77' x 



+ {— + {^Ty + z(z-i) ) v = > 

where (a, a'), (7,7') and (/3,(3') are the exponents at singular points. They satisfy the 
Fuchs relation 

a + a' + 7 + Y + /3 + /3' = l. 

We denote the differences of exponents by 

A = a — a', v = 7 — 7', /i = (3 — j3' . 

Necessary and sufficient conditions for solvability of the identity component of the 
differential Galois group of Q are given by the following theorem due to Kimura 
see also jTTj. 



Theorem 7 The identity component of the differential Galois group of equation (0) is 
solvable if and only if 
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A: at least one of the four numbers A + /i + v, —A + /x + v, A — [i + v, A + fi — v is an 
odd integer, or 

B: the numbers A or —A and fi or —fi and v or —v belong (in an arbitrary order) to 
some of the following fifteen families 



1 


1/2 + 1 


1/2 + to 


arbitrary complex number 




2 


1/2 + 1 


1/3 + 771 


1/3 + q 




3 


2/3 + 1 


1/3 + m 


1/3 + q 


I + m + q even 


4 


1/2 + 1 


1/3 + m 


1/A + q 




5 


2/3 + 1 


1/4 + m 


1/4 + g 


I + m + q even 


6 


1/2 + 1 


1/3 + 771 


1/5 + g 




7 


2/5 + 1 


1/3 + m 


1/3 + g 


I + m + q even 


8 


2/3 + 1 


1/5 + m 


1/5 + g 


I + m + q even 


9 


1/2 + 1 


2/5 + m 


1/5 + g 


I + m + q even 


10 


3/5 + / 


1/3 + m 


1/5 + g 


I + m + q even 


11 


2/5 + / 


2/5 + m 


2/5 + g 


I + m + q even 


12 


2/3 + 1 


1/3 + TO 


1/5 + g 


I + m + q even 


13 


4/5 + / 


1/5 + m 


1/5 + q 


I + m + q even 


U 


1/2 + 1 


2/5 + m 


1/3 + q 


I + m + q even 


15 


3/5 + / 


2/5 + 771 


1/3 + q 


I + m + q even 



Here /, m and q are integers. 
Next we consider the Lame equation in the standard Weierstrass form 

0=[n(n+l)p(t)+B]£, (10) 

where n and B are, in general, complex parameters and p(t) is the elliptic Weierstrass 
function with invariants g2, (?3- In other words, p(t) is a solution of differential equation 

x 2 = f(x), f(x):=4x 3 -g 2 x-g 3 = 4:(x-x 1 )(x-X2)(x-x 3 ). (11) 

We assume that parameters n, B, g 2 and ^3 are such that 

A = g\ - 27g£ ? 0, 

and thus equation f(x) = has three different roots x±, x 2 and x%. All the cases when 
the Lame equation is solvable are listed in the following theorem, see |17j . 

Theorem 8 The Lame equation is solvable only in the following cases 

(i) the Lame and Hermite case (see e.g. IJOlj ) for which n 6 Z and three other 
parameters are arbitrary, 

(ii) the Brioschi-Halphen-Crowford case (see e.g. 14 1[ \40p- In this case 
n + I 6 N and B,g 2 ,gz satisfy an appropriate algebraic equation, 

(Hi) the Baldassarri case ^7J/. Then n + ~ e |Z U -JL U |Z\Z, and there are additional 
algebraic conditions on B,g 2 ,g 3 . 
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Let C(z) denote the set of complex rational functions of z and we consider the second 
order differential equation 



rj" + p(z)rj' + q(z)rj = 0, 



— , p(z),q(z)eC(z) 



Putting 



rj = y exp 



p(s) ds 



we obtain its reduced form 



r(z)y, 



(12) 



For this equation its differential Galois group Q is an algebraic subgroup of SL(2, C). 
The following theorem describes all possible forms of Q and relates them to forms of 
solutions of 1IT2|) . see pilTTj. 

Lemma 1 Let Q be the differential Galois group of equation (112)) . Then one of four 
cases can occur. 

(i) Q is conjugated to a subgroup of the triangular group; in this case equation (|T2*j) has 
a solution of the form y = exp J u>, where u G C(z), 

(ii) Q is conjugated with a subgroup of 



c 
c" 1 



c 6 C* > U 



c 
cr 1 



c e 



in this case equation ()12|) has a solution of the form y = exp J uj, where uj is 
algebraic over C(z) of degree 2, 

(Hi) Q is primitive and finite; in this case all solutions of equation ([12)1 are algebraic, 

(iv) Q = SL(2,C) and equation ()12j) has no Liouvillian solution. 

For a definition of the Liouvillian solution see e.g. (121 • An equation with a Liouvillian 
solution we called integrable. When case (i) in the above lemma occurs, we say that 
the equation is reducible and its solution of the form prescribed for this case is called 
exponential. 

Remark 1 Let us assume that equation ()12|) is Fuchsian, i.e., r(z) has poles at Zi G C, 
i — 1, . . . , K and at zk+\ = oo; all of them are of order not higher than 2. Then at each 
singular point and z = oo we have two (not necessarily different) exponents, see e.g. 
JEBj - One can show, see that an exponential solution which exists when case (i) in 
Lemma[J\ occurs, has the following form 

K 

y = P\[{z-z i y% 
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where is an exponent at singular point; P is a polynomial and, moreover 

K 

deg P = -eoo - ^2 ' 



K 

1=1 

where e M is an exponent at the infinity. 



Remark 2 If equation (|T2*|) has a regular singular point zq with exponents (ei,e2) and 
G\ — &2 Z ; then in a neighbourhood of zq there exist two linearly independent solutions 
of the form 

Vi = (z - z ) ei fi(z), 2 = 1,2, 

where fi{z) are holomorphic at zq. If e\ — e2 G Z ; then one local solution has the 
above form (for the exponent with a larger real part). The second solution can contain 
a logarithmic term, for details see ]H3j- If the logarithmic term appears, then it can be 
shown that only case (i) and case (iv) in LemmaU\can occur, see \2ty . 

3. Particular solution and variational equations 

Without loss of generality, choosing appropriately units of time, mass and length, we 
can put m = g = Iq = 1. Then the Hamiltonian of the generalised spring-pendulum in 
spherical coordinates has the following form 

H = \(p 2 r + ^ + -^K-^ -rcos£ + ^(r-l) 2 -^(r-l) 3 . 



2 V r 2 r 2 sin 2 #; 2 v ; 3 

As we can see ip is a cyclic coordinate and p v is a first integral. Manifold 

J\f = {(r,,0,<p,p ri pe,p v ) e C 6 1 9 = <p = p e = p v = 0} 

is invariant with respect to the flow of Hamilton equations generated by H. Hamiltonian 
equations restricted to M have the form 

r = p r , p r = 1 — k(r — 1) + a(r — l) 2 , 
and can be rewritten as 

r — I — k(r — 1) + a(r — l) 2 . 
Thus the phase curve located on the energy level H\m = E is given by the equation 

E=^+ k -{r-lY- a -{r-lf-r, (13) 

and hence, for the generic values of E, it is an elliptic curve when a ^ (for a = it is 
a sphere). To find its explicit time parametrisation we put 
Q 2a + k 

r = -x+— , 14 

a 2a 

then (|13|) transforms into equation of the form (Jllj) with 

_A; 2 -4a _ k 3 — 6ak — 12a 2 (E + 1) 

92 - — , 93 - ' ^ > 
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For these invariants, x(t) is a non-degenerated Weierstrass function provided that 
(4a - k 2 f + [A; 3 - 6ak - \2a 2 {E + l)] 2 



A 



1728 



7^0. 



But A = only for two exceptional values of energy corresponding to unstable and 
stable equilibria (we assume here that a ^ 0): 

k 3 - 6a{k + 2a) + [k 2 - 4a) 3 / 2 

u — 1 r, o ' 

tf-6 a (k + 2n)-(k 2 -An) 3 / 2 ' 
E K = 



12a 2 

For E B < E < E u we obtain one parameter family T(t,E) of particular solutions 
(r(t), 0, 0,p r (t), 0, 0) expressed in terms of the Weierstrass function and its derivative 

as 



6 k 

r (t) = -p(t;92,g3) + l + — , 

a la 



Pr(t) = -p(t; g2,gs)- 
a 



;i?) 



Particular solutions are single-valued, meromorphic and double periodic with 
periods 2lui and 2a; 2 , and they have one double pole at t — 0. Thus, Riemann surfaces 
F(t, E) are tori with one point removed. 

Using first integral p v , we can reduce the order of VE by two. We choose the zero 
level of this first integral. Let rj = (R, Pr, 0, Pq) denote variations in (r, p r , 6,pg). Then 
the reduced variational equations restricted to the level p v = have the form 
d?/ 



dt 



Lrj, 



where matrix L is given by 



2a(r - 1 




1 

k 















-r 

The normal variational equations read 
1 



(18) 



e 



Pt 



e 



and can be written as 



6 + 2^9 + 4,6 



-r6, 



0. 



(19) 



r{t) r(t) 

Putting $ = Or, and expressing r in terms of the Weierstrass function using ()17|) we 
transform fT9|) to the form 

k 2 - 144p(t) 2 



-$ = 0. 



(20) 



24p(t) + 2k + 4a 

Apart from t = 0, equation (}2T)j) has other singular points which are solutions of the 
equation 

p(t) = d:=~{k + 2a). 
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If d $l {x\, x 2 , X3}, then the above equation has two roots. If d = x^, k = 1,2,3, then 
this equation has one double root. If a = —k, then t = is the only singular point and 
in this case (ffllj) is the Lame equation. 



4. Non-integrability of the classical spring-pendulum 

In this section we investigate the classical spring-pendulum, i.e. we assume that a = 0. 
In this case Hamiltonian equations restricted to manifold M have the form 

r = p r , p r = 1 — k(r — 1), 

and the phase curve corresponding to energy value E is a sphere 

E = 1 — (r — l) 2 -r. 

2 2 V ; 

Making transformation t 1— > z = r(t) we transform NVE (|19|) to a Fuchsian equation 
with rational coefficients and four singular points z — 0, z\ — Z\{E\ z 2 = z 2 (E) and 
23 = cxd, i.e., for generic values of E the transformed NVE is a Heun equation. However, 
changing E we are able to make a confluence of two singular points and for these special 
choices of E the transformed NVE has the form of the Riemann P equation Q. We have 
two possibilities: we can chose E — E\ such that Z\{E\) = z 2 (Ei), or we can take E = E 2 
such that Zi(E 2 ) = 0. In both cases we obtain a Riemann P equation, however these 
two Riemann equations are not equivalent and thus they give two different necessary 
conditions for the integrability. It seems that this fact was not noticed in previous 
investigations. 

Let us assume that k 7^ and put E = —(2k + 1)/ (2k). Then the following change 
of variable 

transforms (fTT)|) to the form 

y " + + T^i) y ' + (- (i + *)!,-!)» + (iw-i) ) y = °' (21) 

where y = y(z) := Q(t(z)). This Riemann P equation has exponents 

a = 0, a' = -1, (3 = 2, = 0, 7 = -7 ' 



VTTk' 

The prescribed choice of the energy corresponds to Ei, i.e., in the generic Heun equation 
two non-zero singular points collapse to one. We prove the following. 

Lemma 2 If k ^ and 

k ^ (m + 2f _1 ' 

where m is a non-negative integer, then equation f!21l) does not possess a Liouvillian 
solution. 
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Proof. Local computation shows that equation (|21j) has logarithms in its formal solutions 
at zero and infinity whenever k ^ 0. Thus, as we know from Remark |21 if the equation 
has a Liouvillian solution, then it must be an exponential one, i.e. we are in case (i) of 
Lemma ^ As the equation is Fuchsian, from Remark^ it follows that such exponential 
solution has the form 



y = z e ° (z-irP(z), 

where is an exponent aA z — i, i — 0,1, and P is a polynomial whose degree m 
satisfies m = —e^ — eo — e\. Moreover, an expansion of an exponential solution of the 
form given above around a singular point does not contain logarithms. However, we 
know that there are formal solutions at z = and z = oo with logarithms. Thus those 
without logarithms corresponds to the maximal exponents, see j3H|- Hence, we must 
put e = a = 0, = (3 = 2, and we may take e± — 7 = 1/a/1 + k. The condition on 
degree of P imposes that 



k — ~ r~r 1 

(m + 2) 2 



with m a non-negative integer. As we excluded such values of k this finishes the proof. 
□ 

Now, for all non-negative integers m, we have (m + 2)~ 2 — 1 < so, as for a physical 
spring we have k > 0, the above lemma shows that equation (|2~T|) has no exponential 
solution (which was the only possible integrable case) and, finally, the NVE is not 
integrable. This ends the proof of Theorem 01 for case a = 0. 

Remark 3 Of course we can prove Lemma using Theorem ^ For equation (|21j) 
differences of exponents are 

A = 1 ' " = 7rW " = 2 ' 

In is easy to notice that case B in the Kimura theorem is impossible. Thus equation (|2*T|) 

is solvable (i.e. the identity component of it differential Galois group is solvable) if and 

only if the condition from case A of the Kimura Theorem is satisfied. The four numbers 

from case A of the Kimura theorem are equal to 

2 2 
\ + H + v = 3+ - -\ + ji + v = 1 + 



VTTk' vT+T 

2 x 2 

A — fi + v = — lH — , A + fi — v = 3 — 



VTTk VTTk 
The condition that at least one of them is an odd integer is equivalent to k = (m+2)~ 2 — 1. 
We gave another proof of Lemma{^ in order to demonstrate a technique which we use 
in the next section. 



To apreciate the relevance of the physical hypothesis k > 0, we prove the following. 
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Lemma 3 If k = (m + 2) 2 — 1 witt m a non-negative integer, then the identity 
component of the differential Galois group of equation (|2*Tjl is Abelian. 

Proof. Proceeding as in the proof of Lemma 121 we conclude that under our assumption 
equation (|2*T|) is solvable if and only if it has a solution of the form y = P/(z — l) m + 2 
with P a polynomial of degree m. Following the method of j23j , we make the change of 
variables y = Y/(z — l) m + 2 in (121 J) and compute the recurrence relation satisfied by the 
coefficients of a power series solution Y — ^2 u n z n at zero. The recurrence is: 

in — m) (n — 2 — m)u n = (n + l)(n + 2)w n+1 . 

The latter always admits a solution such that u_i = u rn+ i — 0, Uq — 1 and u m 7^ 0, 
which proves that for all non-negative integers m the NVE with k = (m + 2)~ 2 — 1 
admits a solution of the form y = P/{z— l) m + 2 with P a polynomial of degree m. 
Thus, the differential Galois group of equation (}2"Tf conjugates to a subgroup of the 
triangular group (case (i) in Lemma P). Moreover, as all exponents are rational, its 
identity component is Abelian. □ 

The above lemma shows that when k < (so, for the negative Young modulus) 
the necessary condition of the Morales- Ramis theory is satisfied for infinite many cases. 
As integrable systems are extremely rare, it is worth checking if, even for non-physical 
values of k excluded in Lemma El the system is integrable or not. 

To answer this question we take E = k/2, and make the following change of variable 

kr(t) 
l ^ Z '- = 2(jfcTl)' 

in equation (|19p. Choosing the prescribed value of energy, we perform a confluence of 
one non-zero singular point with z = in the generic NVE, i.e. this energy corresponds 
to E 2 . The NVE takes the following form 

V " + {h + 2(^T)) y/+ G(l + k){z -I) 2 " 2(1 + ^-1) ) V = °- (22) 

This is exactly the form of NVE which appears in papers [TH IT71 ITS] and the condition 
for its non-integrability is given by ©• Combining the non-integrability conditions for 
equations (}2"Tj) and (}2"2"|) we show the following. 

Theorem 9 The classical spring-pendulum system given by Hamiltonian (QJ) with fcsi 
is integrable only when k = 0. 

Proof. Assume that the system is integrable. Then both NVEs (J21)) and (J22|) are 
integrable, i.e. they possess Liouvillian solutions. Thus, we have 

(m + 2) 2 

for some non-negative integer m, and 
p +p — I 
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But we can rewrite these conditions in the following form 

k = ^—^, s={m-2)\ 

s 

and 

r 1 

k = - , r = -p(p+l). 

1 — r 2 

As we assumed that k ^ 0, both s and r are positive integers. Now, from equality 
1 — s r 
s 1 — r ' 

it follows that r + s — 1, but it is impossible for positive integers r and s. □ 



5. Non-integrability of the generalised spring-pendulum in the case a/0 
and a ^ —k 

NVE given by (J20|) depends on the energy E through the invariants of the Weierstrass 
function, see formula (|15j). The choice of the value of energy is relevant for computation 
and we put 

„ „ 2(3k + 2a)a 2 - 1 

E = hit) \= . 

12a 2 

For this value of the energy we have the following. 

Lemma 4 If a ^ and a ^ —k, then the differential Galois group of the normalized 
NVE (H3J for E = E is equal to SL(2, C). 

Proof. Computation shows that the image of equation (ffllj) under the change of 
variable 1 1— > x — p(t) is the following 

„. , lf'(x) „ . 144x 2 - k 2 - 2a(a + k) . . n , 00 . 

V (x) + --vH-y (x) ; : x „i s y(x) = 0, (23) 

yy! 2f(x) yy! (12x + k + 2a) f(x) KJ K ' 

where 

f(x) = 4 x 3 - g 2 x - g 3 = A(x - Xt)(x - x 2 )(x - x 3 ), 

g 2 and g 3 are given by (|15|) with E = E , and y(x) = $(t(a;)). This equation is Fuchsian 
and it has five singular points: xq = —{2a + k)/12, the three roots xi,x 2 , x 3 of f(x), 
and x 4 = oo. The exponents at the first singularity x are (0, 1), the exponents at the 
roots of f{x) are (0, 1/2), and the exponents at infinity are (—1,3/2). 

If a ^ {0, —k}, then calculation of the formal solutions at xq shows that they 
contain a logarithm. So from Remark |21 we know that the differential Galois group of 
equation is either reducible or it is SL(2, C). 

Let us first assume that the equation is reducible (case (i) of Lemma [TJ, i.e. it has 
an exponential solution. From Remark^ we know that such solution has the form 



P(x)H 



i=0 
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where the e, is an exponent at x = x,, and the degree m of polynomial P(x) satisfies 
_ _ _ 

ifl Coo 2—/i=0 ^i' 

Because the formal solution of valuation at xq has a logarithm, the valuation 
(i.e. the exponent) of y at x must be equal to 1, so e = 1. The exponents at Xi for 
% — 1, 2, 3 of / are of the form n,/2, (rij a non-negative integer), so the relation for the 
degree m is either m = —3/2 — 1 — n/2, or m = 1 — 1 + n/2, for some non-negative 
integer n = ri\ + n 2 + ra 3 . 

If m = —3/2 — 1 — n/2, then m < 0, which is not possible, so we must have 
m = 1 — 1— n/2 = —n/2, which is possible only if n = 0, and hence m = 0. 

So the only possibility is y = (x—xo). Substituting this candidate into equation (J2H|) 
shows that this is not a solution. Hence, the equation is irreducible and, because of the 
logarithms in the local solutions, the only possibility is that the differential Galois group 
is the full SL(2, C), which proves the lemma. □ 

Our main Theorem El in the Introduction now follows, as an immediate consequence 
of the Morales- Ramis theorem, from the considerations in Section 4 (for the case a = 0) 
and from the above lemma. 



6. Non-integrability of the generalised spring-pendulum in the case a = —k 

First we show that for the excluded case a = —k the necessary condition for integrability 
given by the Morales- Ramis theory is fulfilled. 

Lemma 5 For generalised spring-pendulum in case a = —k the identity component of 
differential Galois group of NVE (|20j) is Abelian. 

Proof. For a = —k equation (J2UJ) reads 

* = (6p(t) + ±k\ $, (24) 

so, it has the form of the Lame equation (fTT)|) with n(n + 1) = 6, and B = k/2. For 
the prescribed choice of parameters the invariants g-i and g 3 of the Weierstrass function 
p(t) are following 

_ (k + A)k k 2 (k - 12E - 6) 

92 - — , £3 - ^6 ' 

The discriminant 

_ k 3 [k{k - 12E - 6) 2 - (k + 4) 3 ] 

~ 1728 ' 

is only zero for two exceptional values of E corresponding to two local extrema of the 
potential. Assuming that E is different from these exceptional values, we can apply 
Theorem |H| 

As for equation (|2*4*|) we have n(n + 1) = 6, so n = 2 or n = —3. Since nfZ Lame 
equation is solvable and possesses the Lame-Hermite solutions [El HOI EHj • But for 
a Lame equation with such solutions the differential Galois group is Abelian [T7j . 
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For the excluded energy values, i.e. when E = E s or E = E u (see formula (|16|) ) the 
NVE (after transformation t i— > z := r(t)) has the form of Riemann P equation which 
is solvable. Namely, for both choices of E the case (ii) from Lemma Q occurs. Thus, the 
identity component of the differential Galois group of NVE is Abelian. □ 

Let us notice here that we have at our disposal another family of particular solutions 
corresponding to the following invariant manifold 

TVi = {(r,9,ip,p r ,p 9 ,p v ) G C 6 | ip = p e = p v = 0,9 = n}. 

However, as calculations show, using these particular solutions we do not obtain new 
necessary conditions for the integrability. Because of this, we decide to apply Theorem |3 
Following the decoupling of the first VE into tangential and normal equations (see 
(|18)l). we find that the second variational equations are the following, 

h ~ 12 pr 2 = ( 8 ^ iPe f^ - k (r l} 2 - \ ft) 2 , (25) 
(-12 p + k) 2 



and 



2Ap ■ 2k l6k 3 r lP e A 
+ TToTTTT 2 = 7TT7- 77s 



~12p + k -I2p + k (12 p-k)' 

Wtfpeift 192k 3 p m x k 2 r x 
+ 777 TT3 ~ 777 77i r ^Po,i ~ — 7~2, 



(26) 



[12 p-kf (12 p - k) ' (12 p-kY 

where (r±, 6i,p r ,i-iVe,x) refer to solutions of the first variational system and 
(r2,0 2 ,p r ,2,Pe,2) refer to solutions of the second variational system that we want to 
solve. The equations are now inhomogeneous, with left-hand sides corresponding to the 
(homogeneous) first variational equations, and right-hand sides formed of solutions of 
the first variational equations (which induces coupling). 

These equations look non-linear, at first. However, as explained in ^7], the right- 
hand sides are formed of linear combinations of solutions of the second symmetric powers 
of the first variational system. Hence, the second variational system, together with the 
first, still reduces to a linear differential system and it makes sense to study its differential 
Galois group and its integrability. This fact remains true for variational equations of an 
arbitrary order. 

As the first variational equations are solvable, we could write explicit solutions 
and then solve the second variational equations by variation of constants, but a better 
strategy is to proceed as in [37]: as the first variational equations are Lame equations, 
they have Abelian Galois group if and only if their formal solutions at zero do not 
contain logarithms ^JJ, and it is shown in [SZ] that this remains true for variational 
equations of an arbitrary order. 

This is easily tested in the following way: first we compute formal solutions (as 
a power series) of the first variational equations around zero. Then we plug a generic 
linear combination of these power series in the right hand sides of (j25|) and (|26jh Next 
we apply the method of variation of constants: we thus have to integrate a (known) 
combination of power series and there is a logarithm if and only if this combination of 
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power series has a non-zero residue (i.e. a term of degree —1 in its (Laurent) expansion 
in powers of t). 

Performing this computation we show that the second variational equations are 
integrable. Iterating the process, we computed the solutions of the third, fourth, . . . , 
until the seventh variational equations and found that they are all integrable. We could 
not continue the calculations to higher variational equations for the following reasons. 

The first fact is that the size of the right hand sides of the successive variational 
equations grows rapidly. 

The second fact is that the valuation of the solutions decreases as the order of the 
variational equation grows. For example, the valuations of r2 is —4, the valuation of r^ 
is —5,.., and the valuation of r 7 is —9. To obtain r 6 with an accuracy up to the term 
of degree (to obtain the terms of negative valuation properly, which is all we need 
for integrability by the above remarks), we need to start from an r\ with 27 terms. To 
obtain r 7 , we need to start from an r\ with 30 terms, and so on. The combination of 
these two facts makes the computation intractable for the variational equations of order 
eight. 

The fact that the valuations decrease is no surprise. We know that the restriction 
of the Hamiltonian flow to the invariant manifold Af of Section |3] is an integrable 
system with one degree of freedom. Calculations show that the corresponding solution 
(r, 0, p r , 0) has a valuation at zero that decreases just like the above (and indeed seems 
to govern the lowest valuation in the r»). 

0.8 I 1 1 1 1 1 1 1 




-0.8 1 1 1 1 1 1 1 1 1 

-0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 

Figure 1. Poincare cross section for the generalised spring-pendulum when k = —a = 
4/3 and E = —0.8. The cross section plane with coordinates (8,pe) is fixed at r = 1. 
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Now the fact that the variational equations up to order 7 are integrable might lead 
to a suspicion that the system could be integrable. However, numerical experiments 
clearly indicate chaotic behaviour which contradicts meromorphic integrability. We 
show an example of our numerical experiments in Figure H in this figure we show the 
Poincare cross section for energy E = —0.8 and k = —a = 4/3. On the level H = E we 
chose (r,9,p g ) as coordinates. The cross-section plane was fixed at r = 1. 

The model of the generalised swinging pendulum for a = —k is thus a puzzling 
example of a system that seems (numerically) to be non-integrable but where even a 
deep application of the Morales- Ramis theory is not enough to detect rigorously this 
non-integrability. There is only one reported result concerning application of higher 
variational equations for proving non-integrability. In [37] Morales-Ruiz reports that 
for a certain case of the Henon-Heiles system the identity component of the differential 
Galois group of first and second order VEs is Abelian but for the third order VE it is 
not. We have also several examples of Hamiltonian system for which VE^ have Abelian 
identity component of differential Galois group for k < 3 but non-Abelian for k = 3. 
Thus, as far as we know the generalised spring-pendulum system with a = — k is the 
only example where the application of higher order variational equations is unsuccessful 
in proving non-integrability. 
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